Orthologous Genes Research Articles

Comparative transcriptomics of seed nourishing tissues: uncovering conserved and divergent pathways in seed plants.

The evolutionary and ecological success of spermatophytes is intrinsically linked to the seed habit, which provides a protective environment for the initial development of the new generation. This environment includes an ephemeral nourishing tissue that supports embryo growth. In gymnosperms this tissue originates from the asexual proliferation of the maternal megagametophyte, while in angiosperms it is a product of fertilization, and is called the endosperm. The emergence of these nourishing tissues is of profound evolutionary value, and they are also food staples for most of the world's population. Here, using Orthofinder to infer orthologue genes among newly generated and previously published datasets, we provide a comparative transcriptomic analysis of seed nourishing tissues from species of several angiosperm clades, including those of early diverging lineages, as well as of one gymnosperm. Our results show that, although the structure and composition of seed nourishing tissues has seen significant divergence along evolution, there are signatures that are conserved throughout the phylogeny. Conversely, we identified processes that are specific to species within the clades studied, and thus illustrate their functional divergence. With this, we aimed to provide a foundation for future studies on the evolutionary history of seed nourishing structures, as well as a resource for gene discovery in future functional studies.

Tripartite motif 2b (trim2b) restricts spring viremia of carp virus by degrading viral proteins and negative regulators NLRP12-like receptors.

Tripartite motif (TRIM) proteins are involved in different cellular functions, including regulating virus infection. In teleosts, two orthologous genes of mammalian TRIM2 are identified. However, the functions and molecular mechanisms of piscine TRIM2 remain unclear. Here, we show that trim2b-knockout zebrafish are more susceptible to spring viremia of carp virus (SVCV) infection than wild-type zebrafish. Transcriptomic analysis demonstrates that NOD-like receptor (NLR), but not RIG-I-like receptor (RLR), signaling pathway is significantly enriched in the trim2b-knockout zebrafish. In vitro, overexpression of Trim2b fails to degrade RLRs and those key proteins involved in the RLR signaling pathway but does for negative regulators NLRP12-like proteins. Zebrafish Trim2b degrades NLRP12-like proteins through its NHL_TRIM2_like and IG_FLMN domains in a ubiquitin-proteasome degradation pathway. SVCV-N and SVCV-G proteins are also degraded by NHL_TRIM2_like domains, and the degradation pathway is an autophagy lysosomal pathway. Moreover, zebrafish Trim2b can interfere with the binding between NLRP12-like protein and SVCV viral RNA and can completely block the negative regulation of NLRP12-like protein on SVCV infection. Taken together, our data demonstrate that the mechanism of action of zebrafish trim2b against SVCV infection is through targeting the degradation of host-negative regulators NLRP12-like receptors and viral SVCV-N/SVCV-G genes.IMPORTANCESpring viremia of carp virus (SVCV) is a lethal freshwater pathogen that causes high mortality in cyprinid fish. In the present study, we identified zebrafish trim2b, NLRP12-L1, and NLRP12-L2 as potential pattern recognition receptors (PRRs) for sensing and binding viral RNA. Zebrafish trim2b functions as a positive regulator; however, NLRP12-L1 and NLRP12-L2 function as negative regulators during SVCV infection. Furthermore, we find that zebrafish trim2b decreases host lethality in two manners. First, zebrafish Trim2b promotes protein degradations of negative regulators NLRP12-L1 and NLRP12-L2 by enhancing K48-linked ubiquitination and decreasing K63-linked ubiquitination. Second, zebrafish trim2b targets viral RNAs for degradation. Therefore, this study reveals a special antiviral mechanism in lower vertebrates.

Aspergillus flavus pangenome (AflaPan) uncovers novel aflatoxin and secondary metabolite associated gene clusters

BackgroundAspergillus flavus is an important agricultural and food safety threat due to its production of carcinogenic aflatoxins. It has high level of genetic diversity that is adapted to various environments. Recently, we reported two reference genomes of A. flavus isolates, AF13 (MAT1-2 and highly aflatoxigenic isolate) and NRRL3357 (MAT1-1 and moderate aflatoxin producer). Where, an insertion of 310 kb in AF13 included an aflatoxin producing gene bZIP transcription factor, named atfC. Observations of significant genomic variants between these isolates of contrasting phenotypes prompted an investigation into variation among other agricultural isolates of A. flavus with the goal of discovering novel genes potentially associated with aflatoxin production regulation. Present study was designed with three main objectives: (1) collection of large number of A. flavus isolates from diverse sources including maize plants and field soils; (2) whole genome sequencing of collected isolates and development of a pangenome; and (3) pangenome-wide association study (Pan-GWAS) to identify novel secondary metabolite cluster genes.ResultsPangenome analysis of 346 A. flavus isolates identified a total of 17,855 unique orthologous gene clusters, with mere 41% (7,315) core genes and 59% (10,540) accessory genes indicating accumulation of high genomic diversity during domestication. 5,994 orthologous gene clusters in accessory genome not annotated in either the A. flavus AF13 or NRRL3357 reference genomes. Pan-genome wide association analysis of the genomic variations identified 391 significant associated pan-genes associated with aflatoxin production. Interestingly, most of the significantly associated pan-genes (94%; 369 associations) belonged to accessory genome indicating that genome expansion has resulted in the incorporation of new genes associated with aflatoxin and other secondary metabolites.ConclusionIn summary, this study provides complete pangenome framework for the species of Aspergillus flavus along with associated genes for pathogen survival and aflatoxin production. The large accessory genome indicated large genome diversity in the species A. flavus, however AflaPan is a closed pangenome represents optimum diversity of species A. flavus. Most importantly, the newly identified aflatoxin producing gene clusters will be a new source for seeking aflatoxin mitigation strategies and needs new attention in research.

Chromosome-level Genome Assembly of the Yeast Lodderomyces beijingensis Reveals the Genetic Nature of Metabolic Adaptations and Identifies Subtelomeres as Hotspots for Amplification of Mating Type Loci.

Lodderomyces beijingensis is an ascosporic ascomycetous yeast. In contrast to related species Lodderomyces elongisporus, which is a recently emerging human pathogen, L. beijingensis is associated with insects. To provide an insight into its genetic makeup, we investigated the genome of its type strain, CBS 14171. We demonstrate that this yeast is diploid and describe the high contiguity nuclear genome assembly consisting of eight chromosome-sized contigs with a total size of about 15.1 Mbp. We find that the genome sequence contains multiple copies of the mating type loci and codes for essential components of the mating pheromone response pathway, however, the missing orthologs of several genes involved in the meiotic program raise questions about the mode of sexual reproduction. We also show that L. beijingensis genome codes for the 3-oxoadipate pathway enzymes, which allow the assimilation of protocatechuate. In contrast, the GAL gene cluster underwent a decay resulting in an inability of L. beijingensis to utilize galactose. Moreover, we find that the 56.5 kbp long mitochondrial DNA is structurally similar to known linear mitochondrial genomes terminating on both sides with covalently closed single-stranded hairpins. Finally, we discovered a new double-stranded RNA mycovirus from the Totiviridae family and characterize its genome sequence.

CAM evolution is associated with gene family expansion in an explosive bromeliad radiation.

The subgenus Tillandsia (Bromeliaceae) belongs to one of the fastest radiating clades in the plant kingdom and is characterised by the repeated evolution of Crassulacean acid metabolism (CAM). Despite its complex genetic basis, this water-conserving trait has evolved independently across many plant families and is regarded as a key innovation trait and driver of ecological diversification in Bromeliaceae. By producing high-quality genome assemblies of a Tillandsia species pair displaying divergent photosynthetic phenotypes, and combining genome-wide investigations of synteny, transposable element (TE) dynamics, sequence evolution, gene family evolution and temporal differential expression, we were able to pinpoint the genomic drivers of CAM evolution in Tillandsia. Several large-scale rearrangements associated with karyotype changes between the two genomes and a highly dynamic TE landscape shaped the genomes of Tillandsia. However, our analyses show that rewiring of photosynthetic metabolism is mainly obtained through regulatory evolution rather than coding sequence evolution, as CAM-related genes are differentially expressed across a 24-hour cycle between the two species but are not candidates of positive selection. Gene orthology analyses reveal that CAM-related gene families manifesting differential expression underwent accelerated gene family expansion in the constitutive CAM species, further supporting the view of gene family evolution as a driver of CAM evolution.

A secondary mechanism of action for triazole antifungals in Aspergillus fumigatus mediated by hmg1.

Triazole antifungals function as ergosterol biosynthesis inhibitors and are frontline therapy for invasive fungal infections, such as invasive aspergillosis. The primary mechanism of action of triazoles is through the specific inhibition of a cytochrome P450 14-α-sterol demethylase enzyme, Cyp51A/B, resulting in depletion of cellular ergosterol. Here, we uncover a clinically relevant secondary mechanism of action for triazoles within the ergosterol biosynthesis pathway. We provide evidence that triazole-mediated inhibition of Cyp51A/B activity generates sterol intermediate perturbations that are likely decoded by the sterol sensing functions of HMG-CoA reductase and Insulin-Induced Gene orthologs as increased pathway activity. This, in turn, results in negative feedback regulation of HMG-CoA reductase, the rate-limiting step of sterol biosynthesis. We also provide evidence that HMG-CoA reductase sterol sensing domain mutations previously identified as generating resistance in clinical isolates of Aspergillus fumigatus partially disrupt this triazole-induced feedback. Therefore, our data point to a secondary mechanism of action for the triazoles: induction of HMG-CoA reductase negative feedback for downregulation of ergosterol biosynthesis pathway activity. Abrogation of this feedback through acquired mutations in the HMG-CoA reductase sterol sensing domain diminishes triazole antifungal activity against fungal pathogens and underpins HMG-CoA reductase-mediated resistance.

A Lateral Flow-Recombinase Polymerase Amplification Method for Colletotrichum gloeosporioides Detection

The greater yam (Dioscorea alata), a widely cultivated and nutritious food crop, suffers from widespread yield reduction due to anthracnose caused by Colletotrichum gloeosporioides. Latent infection often occurs before anthracnose phenotypes can be detected, making early prevention difficult and causing significant harm to agricultural production. Through comparative genomic analysis of 60 genomes of 38 species from the Colletotrichum genus, this study identified 17 orthologous gene groups (orthogroups) that were shared by all investigated C. gloeosporioides strains but absent from all other Colletotrichum species. Four of the 17 C. gloeosporioides-specific orthogroups were used as molecular markers for PCR primer designation and C. gloeosporioides detection. All of them can specifically detect C. gloeosporioides out of microbes within and beyond the Colletotrichum genus with different sensitivities. To establish a rapid, portable, and operable anthracnose diagnostic method suitable for field use, specific recombinase polymerase amplification (RPA) primer probe combinations were designed, and a lateral flow (LF)-RPA detection kit for C. gloeosporioides was developed, with the sensitivity reaching the picogram (pg) level. In conclusion, this study identified C. gloeosporioides-specific molecular markers and developed an efficient method for C. gloeosporioides detection, which can be applied to the prevention and control of yam anthracnose as well as anthracnose caused by C. gloeosporioides in other crops. The strategy adopted by this study also serves as a reference for the identification of molecular markers and diagnosis of other plant pathogens.

Intraspecific variation in antibiotic resistance potential within E. coli.

Intraspecific genomic diversity brings the potential for an unreported and diverse reservoir of cryptic antibiotic resistance genes in pathogens, as cryptic resistance can occur without major mutations and horizontal transmission. Here, we predicted the differences in the types of antibiotics and genes that induce cryptic and latent resistance between micro-diverse Escherichia coli strains. For example, we hypothesize that known resistance genes will be the culprit of latent resistance within clinical strains. We used a modified functional metagenomics method to induce expression in eight E. coli strains. We found a total of 66 individual genes conferring phenotypic resistance to 11 out of 16 antibiotics. A total of 14 known antibiotic resistance genes comprised 21% of total identified genes, whereas the majority (52 genes) were unclassified cryptic resistance genes. Between the eight strains, 1.2% of core orthologous genes were positive (conferred resistance in at least one strain). Sixty-four percent of positive orthologous genes conferred resistance to only one strain, demonstrating high intraspecific variability of latent resistance genes. Cryptic resistance genes comprised most resistance genes among laboratory and clinical strains as well as natural, semisynthetic, and synthetic antibiotics. Known antibiotic resistance genes primarily conferred resistance to multiple antibiotics from varying origins and within multiple strains. Hence, it is uncommon for E. coli to develop cross-cryptic resistance to antibiotics from multiple origins or within multiple strains. We have uncovered prospective and previously unknown resistance genes as well as antibiotics that have the potential to trigger latent antibiotic resistance in E. coli strains from varying origins.IMPORTANCEIntraspecific genomic diversity may be a driving force in the emergence of adaptive antibiotic resistance. Adaptive antibiotic resistance enables sensitive bacterial cells to acquire temporary antibiotic resistance, creating an optimal window for the development of permanent mutational resistance. In this study, we investigate cryptic resistance, an adaptive resistance mechanism, and unveil novel (cryptic) antibiotic resistance genes that confer resistance when amplified within eight E. coli strains derived from clinical and laboratory origins. We identify the potential of cryptic resistance genes to confer cross-resistance to antibiotics from varying origins and within multiple strains. We discern antibiotic characteristics that promote latent resistance in multiple strains, considering intraspecific diversity. This study may help detect novel resistance genes and functional genes that could become responsible for cryptic resistance among diverse strains and antibiotics, thus also identifying potential novel antibiotic targets and mechanisms.

OrthoRefine: automated enhancement of prior ortholog identification via synteny

BackgroundIdentifying orthologs continues to be an early and imperative step in genome analysis but remains a challenging problem. While synteny (conservation of gene order) has previously been used independently and in combination with other methods to identify orthologs, applying synteny in ortholog identification has yet to be automated in a user-friendly manner. This desire for automation and ease-of-use led us to develop OrthoRefine, a standalone program that uses synteny to refine ortholog identification.ResultsWe developed OrthoRefine to improve the detection of orthologous genes by implementing a look-around window approach to detect synteny. We tested OrthoRefine in tandem with OrthoFinder, one of the most used software for identification of orthologs in recent years. We evaluated improvements provided by OrthoRefine in several bacterial and a eukaryotic dataset. OrthoRefine efficiently eliminates paralogs from orthologous groups detected by OrthoFinder. Using synteny increased specificity and functional ortholog identification; additionally, analysis of BLAST e-value, phylogenetics, and operon occurrence further supported using synteny for ortholog identification. A comparison of several window sizes suggested that smaller window sizes (eight genes) were generally the most suitable for identifying orthologs via synteny. However, larger windows (30 genes) performed better in datasets containing less closely related genomes. A typical run of OrthoRefine with ~ 10 bacterial genomes can be completed in a few minutes on a regular desktop PC.ConclusionOrthoRefine is a simple-to-use, standalone tool that automates the application of synteny to improve ortholog detection. OrthoRefine is particularly efficient in eliminating paralogs from orthologous groups delineated by standard methods.

Genome-wide identification and evolutionary analysis of the NRAMP gene family in the AC genomes of Brassica species

BackgroundBrassica napus, a hybrid resulting from the crossing of Brassica rapa and Brassica oleracea, is one of the most important oil crops. Despite its significance, B. napus productivity faces substantial challenges due to heavy metal stress, especially in response to cadmium (Cd), which poses a significant threat among heavy metals. Natural resistance-associated macrophage proteins (NRAMPs) play pivotal roles in Cd uptake and transport within plants. However, our understanding of the role of BnNRAMPs in B. napus is limited. Thus, this study aimed to conduct genome-wide identification and bioinformatics analysis of three Brassica species: B. napus, B. rapa, and B. oleracea.ResultsA total of 37 NRAMPs were identified across the three Brassica species and classified into two distinct subfamilies based on evolutionary relationships. Conservative motif analysis revealed that motif 6 and motif 8 might significantly contribute to the differentiation between subfamily I and subfamily II within Brassica species. Evolutionary analyses and chromosome mapping revealed a reduction in the NRAMP gene family during B. napus evolutionary history, resulting in the loss of an orthologous gene derived from BoNRAMP3.2. Cis-acting element analysis suggested potential regulation of the NRAMP gene family by specific plant hormones, such as abscisic acid (ABA) and methyl jasmonate (MeJA). However, gene expression pattern analyses under hormonal or stress treatments indicated limited responsiveness of the NRAMP gene family to these treatments, warranting further experimental validation. Under Cd stress in B. napus, expression pattern analysis of the NRAMP gene family revealed a decrease in the expression levels of most BnNRAMP genes with increasing Cd concentrations. Notably, BnNRAMP5.1/5.2 exhibited a unique response pattern, being stimulated at low Cd concentrations and inhibited at high Cd concentrations, suggesting potential response mechanisms distinct from those of other NRAMP genes.ConclusionsIn summary, this study indicates complex molecular dynamics within the NRAMP gene family under Cd stress, suggesting potential applications in enhancing plant resilience, particularly against Cd. The findings also offer valuable insights for further understanding the functionality and regulatory mechanisms of the NRAMP gene family.

A roadmap of phylogenomic methods for studying polyploid plant genera

AbstractPhylogenetic inference of polyploid species is the first step towards understanding their patterns of diversification. In this paper, we review the challenges and limitations of inferring species relationships of polyploid plants using traditional phylogenetic sequencing approaches, as well as the mischaracterization of the species tree from single or multiple gene trees. We provide a roadmap to infer interspecific relationships among polyploid lineages by comparing and evaluating the application of current phylogenetic, phylogenomic, transcriptomic, and whole‐genome approaches using different sequencing platforms. For polyploid species tree reconstruction, we assess the following criteria: (1) the amount of prior information or tools required to capture the genetic region(s) of interest; (2) the probability of recovering homeologs for polyploid species; and (3) the time efficiency of downstream data analysis. Moreover, we discuss bioinformatic pipelines that can reconstruct networks of polyploid species relationships. In summary, although current phylogenomic approaches have improved our understanding of reticulate species relationships in polyploid‐rich genera, the difficulties of recovering reliable orthologous genes and sorting all homeologous copies for allopolyploids remain a challenge. In the future, assembled long‐read sequencing data will assist the recovery and identification of multiple gene copies, which can be particularly useful for reconstructing the multiple independent origins of polyploids.

Monkeypox virus genomic accordion strategies.

The 2023 monkeypox (mpox) epidemic was caused by a subclade IIb descendant of a monkeypox virus (MPXV) lineage traced back to Nigeria in 1971. Person-to-person transmission appears higher than for clade I or subclade IIa MPXV, possibly caused by genomic changes in subclade IIb MPXV. Key genomic changes could occur in the genome's low-complexity regions (LCRs), which are challenging to sequence and are often dismissed as uninformative. Here, using a combination of highly sensitive techniques, we determine a high-quality MPXV genome sequence of a representative of the current epidemic with LCRs resolved at unprecedented accuracy. This reveals significant variation in short tandem repeats within LCRs. We demonstrate that LCR entropy in the MPXV genome is significantly higher than that of single-nucleotide polymorphisms (SNPs) and that LCRs are not randomly distributed. In silico analyses indicate that expression, translation, stability, or function of MPXV orthologous poxvirus genes (OPGs), including OPG153, OPG204, and OPG208, could be affected in a manner consistent with the established "genomic accordion" evolutionary strategies of orthopoxviruses. We posit that genomic studies focusing on phenotypic MPXV differences should consider LCR variability.

Tracing the evolutionary origin of chordate somites in the hemichordate Ptychodera flava.

Metameric somites are a novel character of chordates with unclear evolutionary origins. In the early branching chordate amphioxus, anterior somites are derived from the paraxial mesodermal cells that bud off the archenteron (i.e., enterocoely) at the end of gastrulation. Development of the anterior somites requires FGF signaling, and distinct somite compartments express orthologs of vertebrate non-axial mesodermal markers. Thus, it has been proposed that the amphioxus anterior somites are homologous to the vertebrate head mesoderm, paraxial mesoderm and lateral plate mesoderm. To trace the evolutionary origin of somites, it is essential to study the chordates' closest sister group, Ambulacraria, which includes hemichordates and echinoderms. The anterior coeloms of hemichordate and sea urchin embryos (respectively called protocoel and coelomic pouches) are also formed by enterocoely and require FGF signals for specification and/or differentiation. In this study, we applied RNA-seq to comprehensively screen for regulatory genes associated with the mesoderm-derived protocoel of the hemichordate Ptychodera flava. We also used a candidate gene approach to identify P. flava orthologs of chordate somite markers. In situ hybridization results showed that many of these candidate genes are expressed in distinct or overlapping regions of the protocoel, which indicates that molecular compartments exist in the hemichordate anterior coelom. Given that the hemichordate protocoel and amphioxus anterior somites share a similar ontogenic process (enterocoely), induction signal (FGF), and characteristic expression of orthologous genes, we propose that these two anterior coeloms are indeed homologous. In the lineage leading to the emergence of chordates, somites likely evolved from enterocoelic, FGF-dependent, and molecularly compartmentalized anterior coeloms of the deuterostome last common ancestor.

The alternative enzymes-bearing tunicates lack multiple widely distributed genes coding for peripheral OXPHOS subunits.

The respiratory chain alternative enzymes (AEs) NDX and AOX from the tunicate Ciona intestinalis (Ascidiacea) have been xenotopically expressed and characterized in human cells in culture and in the model organisms Drosophila melanogaster and mouse, with the purpose of developing bypass therapies to combat mitochondrial diseases in human patients with defective complexes I and III/IV, respectively. The fact that the genes coding for NDX and AOX have been lost from genomes of evolutionarily successful animal groups, such as vertebrates and insects, led us to investigate if the composition of the respiratory chain of Ciona and other tunicates differs significantly from that of humans and Drosophila, to accommodate the natural presence of AEs. We have failed to identify in tunicate genomes fifteen orthologous genes that code for subunits of the respiratory chain complexes; all of these putatively missing subunits are peripheral to complexes I, III and IV in mammals, and many are important for complex-complex interaction in supercomplexes (SCs), such as NDUFA11, UQCR11 and COX7A. Modeling of all respiratory chain subunit polypeptides of Ciona indicates significant structural divergence that is consistent with the lack of these fifteen clear orthologous subunits. We also provide evidence using Ciona AOX expressed in Drosophila that this AE cannot access the coenzyme Q pool reduced by complex I, but it is readily available to oxidize coenzyme Q molecules reduced by glycerophosphate oxidase, a mitochondrial inner membrane-bound dehydrogenase that is not involved in SCs. Altogether, our results suggest that Ciona AEs might have evolved in a mitochondrial inner membrane environment much different from that of mammals and insects, possibly without SCs; this correlates with the preferential functional interaction between these AEs and non-SC dehydrogenases in heterologous mammalian and insect systems. We discuss the implications of these findings for the applicability of Ciona AEs in human bypass therapies and for our understanding of the evolution of animal respiratory chain.

Tandem transcription factors PpNAC1 and PpNAC5 synergistically activate the transcription of the PpPGF to regulate peach softening during fruit ripening.

Peach fruit rapidly soften after harvest, a significant challenge for producers and marketers as it results in rotting fruit and significantly reduces shelf life. In this study, we identified two tandem genes, PpNAC1 and PpNAC5, within the sr (slow ripening) locus. Phylogenetic analysis showed that NAC1 and NAC5 are highly conserved in dicots and that PpNAC1 is the orthologous gene of Non-ripening (NOR) in tomato. PpNAC1 and PpNAC5 were highly expressed in peach fruit, with their transcript levels up-regulated at the onset of ripening. Yeast two-hybrid and bimolecular fluorescence complementation assays showed PpNAC1 interacting with PpNAC5 and this interaction occurs with the tomato and apple orthologues. Transient gene silencing experiments showed that PpNAC1 and PpNAC5 positively regulate peach fruit softening. Yeast one-hybrid and dual luciferase assays and LUC bioluminescence imaging proved that PpNAC1 and PpNAC5 directly bind to the PpPGF promoter and activate its transcription. Co-expression ofPpNAC1 and PpNAC5 showed higher levels of PpPGF activation than expression of PpNAC1 or PpNAC5 alone. In summary, our findings demonstrate that the tandem transcription factors PpNAC1 and PpNAC5 synergistically activate the transcription of PpPGF to regulate fruit softening during peach fruit ripening.

Orthology and Paralogy Relationships at Transcript Level.

Eukaryotic genes undergo a mechanism called alternative processing, resulting in transcriptome diversity by allowing the production of multiple distinct transcripts from a gene. More than half of human genes are affected, and the resulting transcripts are highly conserved among orthologous genes of distinct species. In this work, we present the definition of orthology and paralogy between transcripts of homologous genes, together with an algorithm to compute clusters of conserved orthologous and paralogous transcripts. Gene-level homology relationships are utilized to define various types of homology relationships between transcripts originating from the same ancestral transcript. A Reciprocal Best Hits approach is employed to infer clusters of isoorthologous and recent paralogous transcripts. We applied this method to transcripts from simulated gene families as well as real gene families from the Ensembl-Compara database. The results are consistent with those from previous studies that compared orthologous gene transcripts. Furthermore, our findings provide evidence that searching for conserved transcripts between homologous genes, beyond the scope of orthologous genes, is likely to yield valuable information.

PfbZIP85 Transcription Factor Mediates ω-3 Fatty Acid-Enriched Oil Biosynthesis by Down-Regulating PfLPAT1B Gene Expression in Plant Tissues.

The basic leucine zipper (bZIP) transcription factor (TF) family is one of the biggest TF families identified so far in the plant kingdom, functioning in diverse biological processes including plant growth and development, signal transduction, and stress responses. For Perilla frutescens, a novel oilseed crop abundant in polyunsaturated fatty acids (PUFAs) (especially α-linolenic acid, ALA), the identification and biological functions of bZIP members remain limited. In this study, 101 PfbZIPs were identified in the perilla genome and classified into eleven distinct groups (Groups A, B, C, D, E, F, G, H, I, S, and UC) based on their phylogenetic relationships and gene structures. These PfbZIP genes were distributed unevenly across 18 chromosomes, with 83 pairs of them being segmental duplication genes. Moreover, 78 and 148 pairs of orthologous bZIP genes were detected between perilla and Arabidopsis or sesame, respectively. PfbZIP members belonging to the same subgroup exhibited highly conserved gene structures and functional domains, although significant differences were detected between groups. RNA-seq and RT-qPCR analysis revealed differential expressions of 101 PfbZIP genes during perilla seed development, with several PfbZIPs exhibiting significant correlations with the key oil-related genes. Y1H and GUS activity assays evidenced that PfbZIP85 downregulated the expression of the PfLPAT1B gene by physical interaction with the promoter. PfLPAT1B encodes a lysophosphatidate acyltransferase (LPAT), one of the key enzymes for triacylglycerol (TAG) assembly. Heterogeneous expression of PfbZIP85 significantly reduced the levels of TAG and UFAs (mainly C18:1 and C18:2) but enhanced C18:3 accumulation in both seeds and non-seed tissues in the transgenic tobacco lines. Furthermore, these transgenic tobacco plants showed no significantly adverse phenotype for other agronomic traits such as plant growth, thousand seed weight, and seed germination rate. Collectively, these findings offer valuable perspectives for understanding the functions of PfbZIPs in perilla, particularly in lipid metabolism, showing PfbZIP85 as a suitable target in plant genetic improvement for high-value vegetable oil production.

Dissecting the causal polymorphism of the Lr67res multipathogen resistance gene.

Partial resistance to multiple biotrophic fungal pathogens in wheat (Triticum aestivum L.) is conferred by a variant of the Lr67 gene, which encodes a hexose-proton symporter. Two mutations (G144R, V387L) differentiate the resistant and susceptible protein variants (Lr67res and Lr67sus). Lr67res lacks sugar transport capability and was associated with anion transporter-like properties when expressed in Xenopus laevis oocytes. Here, we extended this functional characterization to include yeast and in planta studies. The Lr67res allele, but not Lr67sus, induced sensitivity to ions in yeast (including NaCl, LiCl, KI), which is consistent with our previous observations that Lr67res expression in oocytes induces novel ion fluxes. We demonstrate that another naturally occurring single amino acid variant in wheat, containing only the Lr67G144R mutation, confers rust resistance. Transgenic barley plants expressing the orthologous HvSTP13 gene carrying the G144R and V387L mutations were also more resistant to Puccinia hordei infection. NaCl treatment of pot-grown adult wheat plants with the Lr67res allele induced leaf tip necrosis and partial leaf rust resistance. An Lr67res-like function can be introduced into orthologous plant hexose transporters via single amino acid mutation, highlighting the strong possibility of generating disease resistance in other crops, especially with gene editing.

Meta-Quantitative Trait Loci Analysis and Candidate Gene Mining for Drought Tolerance-Associated Traits in Maize (Zea mays L.).

Drought is one of the major abiotic stresses with a severe negative impact on maize production globally. Understanding the genetic architecture of drought tolerance in maize is a crucial step towards the breeding of drought-tolerant varieties and a targeted exploitation of genetic resources. In this study, 511 quantitative trait loci (QTL) related to grain yield components, flowering time, and plant morphology under drought conditions, as well as drought tolerance index were collected from 27 published studies and then projected on the IBM2 2008 Neighbors reference map for meta-analysis. In total, 83 meta-QTL (MQTL) associated with drought tolerance in maize were identified, of which 20 were determined as core MQTL. The average confidence interval of MQTL was strongly reduced compared to that of the previously published QTL. Nearly half of the MQTL were confirmed by co-localized marker-trait associations from genome-wide association studies. Based on the alignment of rice proteins related to drought tolerance, 63 orthologous genes were identified near the maize MQTL. Furthermore, 583 candidate genes were identified within the 20 core MQTL regions and maize-rice homologous genes. Based on KEGG analysis of candidate genes, plant hormone signaling pathways were found to be significantly enriched. The signaling pathways can have direct or indirect effects on drought tolerance and also interact with other pathways. In conclusion, this study provides novel insights into the genetic and molecular mechanisms of drought tolerance in maize towards a more targeted improvement of this important trait in breeding.

A de novo genome assembly of Solanum bulbocastanum Dun., a Mexican diploid species reproductively isolated from the A-genome species, including cultivated potatoes.

Potato and its wild relatives are distributed mainly in the Mexican highlands and central Andes of South America. The South American A-genome species, including cultivated potatoes, are reproductively isolated from Mexican diploid species. Whole-genome sequencing has disclosed genome structure and similarity, mostly in cultivated potatoes and their closely related species. In this study, we generated a chromosome-scale assembly of the genome of a Mexican diploid species, Solanum bulbocastanum Dun., using PacBio long-read sequencing, Optical mapping, and Hi-C scaffolding technologies. The final sequence assembly consisted of 737.9 Mb, among which 647.0 Mb were anchored to the 12 chromosomes. Compared with chromosome-scale assemblies of S. lycopersicum (tomato), S. etuberosum (non-tuber-bearing species with E genome), S. verrucosum, S. chacoense, S. multidissectum, and S. phureja (all four are A-genome species), the S. bulbocastnum genome was the shortest. It contained fewer transposable elements (56.2%) than A-genome species. A cluster analysis was performed based on pairwise ratios of syntenic regions among the seven chromosome-scale assemblies, showing that the A-genome species were first clustered as a distinct group. Then, this group was clustered with S. bulbocastanum. Sequence similarity in 1,624 single-copy orthologous gene groups among 36 Solanum species and clones separated S. bulbocastanum as a specific group, including other Mexican diploid species, from the A-genome species. Therefore, the S. bulbocastanum genome differs in genome structure and gene sequences from the A-genome species. These findings provide important insights into understanding and utilizing the genetic diversity of S. bulbocastanum and the other Mexican diploid species in potato breeding.